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Abstract 

The purpose of this paper is to study the effects of dark energy on 
dynamics of the collapsing fluid within the framework of metric f(R) 
gravity. The fluid distribution is assumed to be locally anisotropic and 
undergoing dissipation in the form of heat flow, null radiations and 
shear viscosity. For this purpose, we take general spherical symmetric 
spacetime. Dynamical equations are obtained and also some special 
solutions are found by considering shearing expansionfree evolution of 
the fluid. It is found that dark energy affects the mass of the col- 
lapsing matter and rate of collapse but does not affect the hydrostatic 
equilibrium. 

Keywords: f{R) theory; Dissipative anisotropic fluid; Dynamical equa- 
tions. 
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1 Introduction 

Dark energy (DE) and gravitational collapse are the two noteworthy issues 
of cosmology and gravitational physics. Recent observational data [T]-[TT] 
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indicate that our universe is expanding. This acceleration is explained in 
terms of DE, which may be explained in modified gravity models. On the 
other hand, gravitational collapse is the basic process driving evolution within 
galaxies, assembling giant molecular clouds and producing stars. 

When the Einstein-Hilbert (EH) gravitational action in General Relativ- 
ity (GR), 
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Seh =2k,J dTxy/^gR, (1.1) 
is re-written in the modified form as follows 

Smodif = y k \ d 4 xyf^f(R), (1.2) 

the presence of f(R) function may be understood as the introduction of an 
effective fluid which is not restricted to hold the usual energy conditions. An 
important feature of this theory is that the modified field equations can be 
written in the form of Einstein tensor which makes it easy to compare with 
GR. This is done by taking all the higher order corrections to the curvature 
on the right hand side of the field equations and defining it as a "Dark source" 
term or " curvature fluid" . In the following, the field equations in the metric 
approach are obtained as 

F{R)R a p~f{R)g a p-V a V p F{R)+g aP nF{R) = KT a p, {a, $ = 0, 1, 2, 3), 

(1.3) 

where F(R) = df(R)/dR. When we re-write this equation in the above 
mentioned form, it follows that 



<7«* = £TO + ?5), (1-4) 
f(R) + RF(R) 



where 



rpD _ 1 



-g a p + V a VpF{R)-g a pUF{R) 



(1.5) 



In this way, DE can be thought of as having the geometrical origin rather 
than some additional scalar fields which are added by hand to the matter 
part. Therefore, f(R) theory of gravity may be used to explain the present 
accelerating expansion of the universe. 

This theory has many applications in cosmology and gravity such as in- 
flation, local gravity constraints, cosmological perturbations and spherically 
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symmetric solutions in weak and strong gravitational backgrounds. In the 
last few years, a considerable amount of analysis and theoretical observations 
have been made in order to compare the preliminary successes of f(R) mod- 
els with the great achievements of GR. For example, Capozziello et al. [T2] 
analyzed the relation between spherical symmetry and the weak field limit 
of f(R) theory and compared the results with GR. In strong gravitational 
background such as neutron star and white dwarfs, one needs to take into 
account the backcreation of gravitational potentials to the field equations. 
The structure of the relativistic stars in f(R) theory has been discussed by 
many authors [T3]-[TT]. 

Cai et al. [18J derived the generalized Misner-Sharp energy in f(R) grav- 
ity for spherically symmetric spacetime. They found that unlike GR, the 
existence of the generalized Misner-Sharp energy depends on a constraint 
condition. Erickcek et al. [19] found unique exterior solution for a stellar 
object by matching it with interior solution in the presence of matter source. 
Kainulainen et al. [20] studied the interior spacetime of stars in Palatini 
f(R) gravity, de la Cruz-Dombriz et al. [21J discussed the problem of find- 
ing static spherically symmetric black hole solutions in f(R) theory. They 
explored several aspects of constant curvature solutions and thermodynam- 
ical properties. In a recent paper [22], we have investigated perfect fluid 
gravitational collapse in this theory and found that constant scalar curva- 
ture term f(Ro) acts as a source of repulsive force and thus slows down the 
collapse of matter. 

In GR, Oppenheimer and Snyder |23j innovated the first mathematical 
model for the description of gravitational collapse of stars. After that many 
approaches are adopted for the physical description of the fluid in order to 
form self-gravitating objects. During fluid evolution, self-gravitating objects 
may pass through phases of intense dynamical activities for which quasi-static 
approximation is not reliable. For instance, the collapse of very massive stars 
[21], the quick collapse phase yielding neutron star formation [25] and the 
peculiar stars. The peculiar stars are very dense, strongly magnetic and are 
created when massive stars die by collapse. Many of the cooler chemically 
peculiar stars are the result of the mixing of nuclear fusion products from 
the interior of the star to its surface. 

Misner and Sharp [26] discussed the gravitational collapse by taking 
spherically symmetric ideal fluid. They provided a full account of the dy- 
namical equations governing the adiabatic relativistic collapse. Dissipative 
process plays dominant role in the formation and evolution of stars. Vaidya 
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[2T] introduced the idea of outgoing radiations in collapse giving non-vacuum 
exterior outside the stars. It was physically a quite reasonable assumption as 
radiation is a confirmation that dissipative processes are occurring, causing 
loss of thermal energy of the system which is an effective way of decreasing 
internal pressure. Cai and Wang |28j studied the formation of black holes 
in the background of DE. Herrera et al. [25] investigated the dynamics of 
gravitational collapse which undergoes dissipation in the form of heat flow 
and radiation. The same authors [30] also provided detailed discussion on 
the physical meaning of expansionfree fluid evolution. Di Prisco et al. [31] 
explored gravitational collapse by adding charge and dissipation in the form 
of shear viscosity. The dynamical and transport equations were coupled to 
observe the effects of dissipation over collapsing process. Recently, Sharif et 
al. [32j-[39j explored different aspects of gravitational collapse by using all 
the three types of symmetry. 

In this paper, we discuss how DE generated by curvature fluid affects the 
dynamics of dissipative gravitational collapse. The format of the paper is 
as follows. In next section 2, we present spacetimes and energy- momentum 
tensor for dissipative fluid. Section 3 is devoted to formulate the modified 
field equations and dynamical equations in f(R) gravity. In section 4, special 
solutions are discussed. The last section 5 concludes the main results of the 
paper. 

2 Spacetimes and Collapsing Matter 

We take spherical symmetry about an origin O which is divided into two 
regions, interior and exterior by 3D hypersurface £ centered at O. The 
interior spacetime to £ can be represented by the line element 



ds 2 = A 2 (t, r)dt 2 - B 2 (t, r)dr 2 - C 2 (t, r){dB 2 + sin 2 ddtf). (2.1) 



For the exterior spacetime to E, we take the Vaidya spacetime given by the 
line element 



where m[y) represents the total mass and v is the retarded time. In the 
interior region, we assume a distribution of anisotropic collapsing fluid which 





(2.2) 



r 
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undergoes a dissipation in the form of heat flow, null radiations and shearing 
viscosity. The energy-momentum tensor with such properties is given by 

T af 3 = {p+P±)u a up-p ± g a p+(p r -p 1 _)x a Xi3+u a q/3+q a up+elJ l3 -2ri<j a p. (2.3) 

Here, we have p as the energy density, p± the tangential pressure, p r the 
radial pressure, q a the heat flux, rj the coefficient of shear viscosity, u a the 
four-velocity of the fluid, Xa the unit four-vector along the radial direction, l a 
a radial null four-vector and e the energy density of the null fluid describing 
dissipation in the free streaming approximation. These quantities satisfy the 
relations 

U U a 1, \ Xa 1) I U a 1, 

u a q a = 0, X a u a = 0, l a l a = (2.4) 
which are obtained from the following definitions in co-moving coordinates 
u a = A~ 1 5^ x a = B- 1 5^, q a = qB- 1 5?, l a = A~ l 5* + B~ l 5« . (2.5) 
Here q is a function of t and r. The shear tensor is defined by 

°al3 = «(a;/3) ~ «(a«/9) ~ ^©(fl'a/J ~ U a Up), (2.6) 

where the acceleration a a and the expansion are given by 

a a = UatfU 13 , 6 = u* a . (2.7) 

The bulk viscosity does not appear explicitly as it has been absorbed in the 
form of radial and tangential pressures of the collapsing fluid. From Eqs.( l2.5j ) 
and (I2.6p . the non-zero components of the shear tensor are 

2 1 
a n = --B 2 a, cr 2 2 = -C 2 a, a 33 = cr 22 sin 2 9. (2.8) 

The shear scalar a is given by 

1 ( B C\ 

Using Eqs.(ES]) and (1277)1 . it follows that 

fl i = -^ a2 = aQaa= (^)' = il^ + 2 ^l' ( 2 - 10 ) 
where dot and prime represent derivative with respect to t and r respectively. 
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3 The Field Equations and the Dynamical 
Equations in f(R) Gravity 

The field equations (11 .3p for the interior metric take the following form 

AA" _B AB AA'B' _ 2C 2AC 2AA'C _ A 2 f(R) 
~W ~B + AB B~ 3 ~C + ~AC + B 2 C 2 V 
A 2 F" F ,—B 2C\ F'A 2 -B' 2C 



B 2 F F^B + ^ ] - FW<^ + = 8 * ( " + e)A • (31) 

,C' CA' BC\ F' A'F BF' , . . 

" 2( C - CA ~ BC > " F + AF + BF = + ' ] *"< ^ 

-A" BB ABB A'B' 2C" 2CBB 2B'C 



A A 2 A 3 AB C A 2 C BC 

B 2 f(R) B 2 F FB 2 ,A 2C , F'A' 2C\ 

H ^^—^ 1 f 1 ) H f 1 ) 

2 F A 2 F FA 2 A C F A C 

= 8n(p r + e + -w)B 2 , (3.3) 

C AC A'C C" CB B'C 1 

+ 



CA 2 CA 3 B 2 AC B 2 C A 2 BC B 3 C C 2 
C 2 C' 2 C 2 f(R) F F" F A B C, 



A 2 C 2 B 2 C 2 2 F A 2 B 2 A 2 A B C 

F'A' B' C\ , 2 . . , 

~ B + C ] = S ^ " a"'" (3 ' 4) 

In the most general case, we have nine variables with four equations. One 
cannot solve this system of equations unless some assumptions are imposed. 
An expansionfree motion of the system, i.e., = 0, will be used in section 4 
which may lead to some interesting results. 

In order to develop dynamical equations that help to study the properties 
of collapsing process, we shall use Misner and Sharp [26] formalism. The mass 
function is defined by 

C CI C 2 C' 2 \ 

M = g (1 + = - I 1 + - 2 - — 2 J . (3.5) 

The proper time and radial derivatives are given by 
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where C is the areal radius of a spherical surface inside the boundary. The 
velocity of the collapsing fluid is defined by the proper time derivative of C, 
i.e., 



U = D T C=j, 



(3.7) 



which is always negative. Using this expression, Eq.f l3.5p implies that 

2M /2 



E 



a_ 

B 



[1 + u 2 + 



c 



(3.8) 



When we make use of Eqs.( P , (g3ED and (JHJ| in Eq.flHZ2), we obtain 



E[~D c (e - a) - 



47T 



(q + e)- + D T F' - ^-D C AD T F - ~^D C B 



. (3.9) 



The rate of change of mass in Eq. (13.51) with respect to proper time, with 
the use of Eqs. fl3.ip - fl3.4l) . is given by 



D T M = — 



~F 



- U 



2C 2 



-4vr{(p + 2e + p r - 2p L + -r)(j)U + E(q + e)} 

o 



c 2 



+ 



ED A D T F f(R)\ 
2C 2 

F' 

)D C A + 




2C 



C 2C" ° 2SC" 
7 [D T F' + ED C FD T B] 




F F F f D T A 
h — + — (— ^ 

C 2 C l 2C 



(3.10) 



This represents variation of total energy inside a collapsing surface of radius 
C. The first two terms inside the square brackets have negative signs which 
show that the total energy is being dissipated during collapse. However, in 
the case of collapse U < 0, the first term (p + 2e + p r — 2p± + f^o") increases 
the energy density through the rate of work being done by the effective 
anisotropic pressure and radiation density of the null fluid. Here we may use 
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equation of state to change energy density into pressure. The second term 
E(q + e) has negative sign which describes that energy is leaving the system 
due to heat flux and radiations. All other terms show the contribution of 
the DE in the form of function f(R) and its derivatives. We know that DE 
exerts a repulsive force on its surrounding thus we may conclude from the 
above expression that DE reduces the mass of the collapsing matter due to 
its negative pressure. It is mentioned here that in GR, only the first two 
terms excluding energy density p and tangential pressure appear. 
Similarly, we can calculate 

C 2 
2F , 

^ - d t fd c a\ - u [ed c f f^- 

EFDnA 1 „ „ F"E 2 ED T F' 



D C M 



8n{p + 2e + Pr + ^qa + -{q + e)} 



+ U 2 



ED C A 
C 



+ 



'Cj 

CO 



+ E 2 D C FD C B) 



FE 2 



C 2 



C' 2 
F 

+ C 2 + 



o 

2FD T U 
C 



^(D T FD T B 



(3.11) 



This equation describes how different quantities influence the mass between 
neighboring surfaces of radius C in the fluid distribution. The first two terms 
and their description in the above expression are almost the same as in GR 
[30] except for the factor {p r + ^rj). The appearance of this factor is due to 
the complicated field equations. The remaining terms represent contribution 
of DE due to curvature fluid. Taking integral of Eq. (13. lip over C, we have 



M 



C^2 

2F 



4 U 
8tt{ P + 2e + p r + -rjcr + -(q + e)} 



+ U 



+ 



c 

EFD C A 
CO 



- D T FD C A 



U < ED C F 



D T B ED C A 



O 



C 



Dj'j'F 



+ E 2 D C FD C B) + 



FE 2 

-Q2 



F"E 2 
" c ,2 ~ 

F 

+ C 2 + 



- - -(D T FD T B 

2FD T U 



C 



o 

dC. 



(3.12) 



The dynamical equations can be obtained from the contracted Bianchi 
identities. Consider the following two equations 



T a0 u 
~0 



0. 



T a0 Y 

;/3 A.a 







(3.13) 
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which yield 



1 

A 
1 

+ B 



A 



4 B 2 C 

(p + e) ■ + (p + 2e + p r + -rjF) - + 2(p + e + P± - - V a) - 



(q + e)' + 2(q + e) 



(Acy 

AC 



(3.14) 



(g + e )+2(g + e)(| + ^; 



1 

B 



4 A! C 
+ {p + p r + 2e + -w)- + 2(p r - P± + e + 2r/a) ^ 

Using Eqs.flZSD, fl2TTU|) . f l3T6|) and (|51ty . it follows that 



0. 



(3.15) 



D T (p + e) + i(3p + 4e + p r + 2p ± )6 + -(e + p r - p ± - 2pa)a 

+SD c (g + e) + 2(q + e)(a + ^) = 0, (3.16) 

2 4 
D T (q + e) + - (q + e)(26 + a) + ££>c(Pr + e + t^) + (p + Pr 

4 £ 
+2e + -pa)a + 2(p r -p± + e + 2r]a) — = 0. (3.17) 

The acceleration DtU of the collapsing matter inside the hypersurface is 
obtained by using Eqs. (l3T4"|l - (l3TB|l and (JUS) 

2M 8vrC. 2 . ^ ^ ^ UED T B 
D T U = — + —(p ± -- W ) + Ea + ED c E 

- ^^--j{-D TT F + ED c (^)-aD T F 



+ e^ + L )Dc f } . 



(3.18) 
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Substituting a from Eq. (13.18j ) into (J3TTTJ) , it follows that 



4 4 
(p + Pr + 2e + -w)D T U = (p + p r + 2e + -rja) 



8nC . 2 , 



2M Cf(R) C 
-^T--^p 1 + pD TT F-aD T F 



4 2 
£>c(Pr + e + -7/a) + — (p r 



— E 

(UD C A 1 \ CD C F 



s , 4 s [UD C A 1 \ 

-p ± + e + 2r/a) + (p + p r + 2 e + -r ?( T) ( — ^— + - J 



-E 



2U 4 
Dt(q + e) + 2(g + e)(— + a) + (p + p r + 2e + -r/a) 



^ ^ UD T B C ^ F\ 
yi-DcE + -^ + -D c (- 



(3.19) 



This shows the role of different forces on the collapsing process. The term 
within the brackets on the left hand side stands for "effective" inertial mass 
and the remaining term is acceleration. The first term on the right hand side 
represents gravitational force (the passive gravitational mass by equivalence 
mass). The term within the first square brackets shows how dissipative terms 
and DE affect the passive gravitational mass. The first two terms in the 
second square brackets are gradient of the effective pressure and effect of 
local anisotropy of pressure with negative sign which increases the rate of 
collapse. While the other term shows the contribution of DE collective with 
effective pressure. Here the presence of U (with negative sign) indicates that 
this term slows down the rate of collapse. The last square brackets depicts 
the combine role of each type of matter component. 



4 Some Special Solutions 

4.1 Shearing Expansionfree Dissipative Fluid 

Now we use the condition of expansionfree motion, i.e., G = 0, to find some 
solutions. Using this condition, Eq. fl2.10p yields 



B 
B 



-2°-. 

c 



(4.1) 
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The physical meaning of this condition is discussed with detail in [3U]. On 
integration, we get 

B= 9 ~P. (4.2) 
where gi(r) is an arbitrary function. Substituting Eq.(4.1) in (13. 21) . we get 

,C> A' C\ ,F' A'F^C F' . ,ABC , . 

2< -C-A- 2 C > + ~F - AF^C + 2 T = 81(9 + (43) 
which, by integrating, gives 

A = F2dC \ - 1 (Mff+O^r+^-^Jg* ( 4 4) 
r(t) 



Thus the interior metric becomes 



2 



-{%fdr 2 - C 2 (t,r)(d6 2 + sin 2 9d<j) 2 ). (4.5) 

The above metric represents spherically symmetric anisotropic fluid which is 
going shearing expansionfree evolution. Assuming the condition of constant 
scalar curvature (R = R c ), according to which F(R C ) = constant, the metric 
for nondissipative (q = e = 0) case is reduced to 

( 77-2/^2 \ 2 1 
dt 2 -— dr 2 -C 2 {t, r){d0 2 + sin 2 8d<f) 2 ). (4.6) 
T{t) J C 4 

Here we take g\(r) = 1. It is mentioned here that f(R) theory, in the 
case of constant scalar curvature, exhibits behavior just like solutions with 
cosmological constant in GR. This is one of the reason why the DE issue 
can be addressed using this theory. Moreover, in this case, Birkhoff theorem 
holds, i.e., stationary solutions are also static p2] which does not hold for 
every f(R). 



4.2 Shearing Expansionfree Nondissipative Perfect Fluid 

Here we discuss shearing nondissipative expansionfree perfect fluid with con- 
stant scalar curvature. For perfect fluid n = and hence the metric reduces 
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to Eq. (j4.6p . Using values of A and B from Eq. (I4.6|) in the field equations 
Eqs.(JXI])-(IS3D along with Eq.flUSP, it follows that 

(\<r 2 n4n" 
() ' + ^J_ + C 3 (2C" 



F*C G C 



^ : ^_i )+c ^ + 4C "' 



F C 4 CC 5V C t' y CC 



8C ' C ') + 10C' 2 C 2 

c 


1 f(Rc) 

2 F c 


8vrp 




D c cr o 
+ 

3 C 


= o, 




C'C IOC' 2 , 

cc c 2 ' 


, i f(Rc) 

2 F c 


8lip r 
1 c 


3C' 2 C" 1 


, i f(Rc) 


87rpj_ 


c 2 c'^c 2 


2 F c 


F c 



F C 4 CC 5 C r 7 CC 

In this case, the Bianchi identities are reduced to the same expression as in 
GR, i.e., 

p + 2 ( p± -p r )-=0, (4.7) 

P' r + {P + Pr)^ + 2(p + 2p r -P±)^ = 0. (4.8) 

We would like to mention here that for isotropic fluid, i.e., p r = pj_, Eq. (14.71) 
implies that energy density p depends only on r. 



4.3 Shearing Expansionfree Dust 

In this case, we have p r = = p±_ and hence Eq. fl4.8p becomes 



whose integration gives 



^ + 2-^0 (4.9) 



6=^fi, C = «&, (4.10) 



where g2{t) and g${r) are arbitrary functions of t and r respectively. Further, 
integrating Eq. (14.101) with respect to the corresponding arguments of the 
arbitrary functions, we get 

C 3 

' 92{t)dt+ I g 3 (r)dr (4.11) 



3 
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which can be written as 

C 3 = m + X(r), (4.12) 

where 

ij(t) = 3 J g 2 {t)dt, X {r) = 3 J g 3 (r)dr, (4.13) 

In view of Eg. (14. 10 p and taking r(t) = Q2{t), the metric in Eq. fl4.6p is reduced 
to 

ds 2 = F*dt 2 - -^dr 2 - C 2 (t, r)(d6 2 + sin 2 6d(j) 2 ) (4.14) 
while Eq.( )3.7p becomes 

U=£ (4.15) 

c 

When we use the standard GR limit, i.e., F c — 1, Eqs.f l4.l4j) and (I4.15P are 
reduced to the corresponding GR results. 



5 Summary 

This paper investigates the gravitational collapse of a spherically symmetric 
star in f(R) theory of gravity. The star is made up of viscous anisotropic 
fluid distribution which is dissipating energy in the form of heat flow, null 
radiations and shearing viscosity. The objective of this work is to explore the 
effects of DE which is generated by modifying EH action. We have discussed 
the consequences in view of existing GR results. 

It is concluded that the contribution of DE terms decreases the mass of 
the collapsing fluid with the passage of time and hence prevents the fluid to 
collapse. In dynamical equations, such curvature terms appear oftentimes to 
affect the passive gravitational mass and rate of collapse. 

We have found some solutions by assuming that fluid has no expansion, 
i.e., = 0. Herrera et al. [30] has made a comprehensive discussion on 
the physical meaning of expansionfree evolution of the fluid. According to 
them such an assumption causes the formation of a vacuum cavity inside 
the collapsing fluid. Using this condition, we have obtained general metric 
for dissipative fluid which is further reduced to nondissipative case. For 
nondissipative solution, a locally anisotropic perfect fluid solution is obtained 
with the assumption of constant scalar curvature. This is further reduced to 
the dust case. 
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It is noted here that implication of hydrostatic equilibrium limit yields 
the same expression as in GR [30J, i.e., 



D C Pr + 2 



Or ~ Rl) 

c 



C{C - 2M) 



(P + Pr) 



(M + A7r Pr C 3 ). 



(5.1) 



The reason is that f(R) term does not affect contracted Bianchi identities 
which are used to obtain the above expression along with the use of Eq. (14.7ft 
giving time independent p in an isotropic pressure case. This change in 
energy density can be interpreted as the rate of work being done by the force 
of locally anisotropic pressure. It would be worthwhile to investigate these 
issues using other symmetries in f(R) theory for the complete understanding 
of gravitational collapse. 
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